function  Hv  = H(a,b)
%Возвращает вектор фокусировки для нашей антенной сиситемы
%Входные аргументы b-азимут и a-угол места.


% Nant=4; %число антенн
%координаты антенн
r{1}=[1/4,-1/4,0];%координаты первой антенны
r{2}=[-1/4,-1/4,0];
r{3}=[-1/4,1/4,0];
r{4}=[-1/4,-1/4,0];%считаем сразу в лямбдах

Nant=length(r);



[kx,ky,kz] = sph2cart(b, a, 1);
k=[kx;ky;kz];
phi = nan(Nant,1);
for i=1:Nant
    
    phi(i)=2*pi*r{i}*k;
    
end

Hv = exp(1i*phi);

end

